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Abstract 

We consider general curvature-invariant modifications of the Einstein-Hilbert action that become 
important only in regions of extremely low space-time curvature. We investigate the far future 
evolution of the universe in such models, examining the possibilities for cosmic acceleration and 
other ultimate destinies. The models generically possess de Sitter space as an unstable solution 
and exhibit an interesting set of attractor solutions which, in some cases, provide alternatives to 
dark energy models. 
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I. INTRODUCTION 



The acceleration of the universe presents one of the greatest problems in theoreti- 
cal physics today. The increasingly accurate observations of type la supernovae light- 
curves^ coupled with exquisite measurements of CMB anisotropics and large scale structure 
data have forced this issue to the forefront of those facing particle physicists, 

cosmologists and gravitational physicists alike. 

This problem has been attacked head on, but no compelling, well-developed and well- 
motivated solutions have yet emerged. While much work has focused on the search for new 
matter sources that yield accelerating solutions to general relativity, more recently some 
authors have turned to the complementary approach of examining whether new gravitational 
physics might be responsible for cosmic acceleration. 

There have been a number of different attempts to modify 

gravity to yield accelerating cosmologies at late times. The path we are concerned with 
in this paper is the direct addition of higher order curvature invariants to the Einstein- 
Hilbert action. The first example of this was provided by the model of Carroll, Duvvuri, 
Trodden, and Turner (CDTT) HI (see also M)- For subsequent work on various aspects and 



extensions of this model, see 




• subsequent worK on various aspects anc 

3,0,0,0,5, 22, 3 0,0 HQ 



In particular, the simplest model has been shown to conflict with solar system tests of 
gravity 

000- 

Our approach is purely phenomenological. The evidence for cosmic 
acceleration is very sound. In pure Einstein gravity, as matter dilutes away in the expanding 
universe, the expansion rate inevitably slows. Our question is can this be avoided within the 
gravitational sector of the theory, as opposed to adding new energy sources. One way is by 
adding a cosmological constant. In this paper we explore a wider class of modifications that 
share the feature of late-time accelerating behavior, thus fitting cosmological observations. 
Consider a simple correction to the Einstein-Hilbert action, 

where // is a new parameter with units of [mass] and Cm is the Lagrangian density for matter. 
This action gives rise to fourth-order equations of motion. In the case of an action depending 



d 4 x y/—g C 
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purely on the Ricci scalar (and on its derivatives) it is possible to transform from the frame 
used in (jTj), which we call the matter frame, to an Einstein frame, in which the gravitational 
Lagrangian takes the Einstein-Hilbert form and the additional degrees of freedom (H and 
H) are represented by a scalar field <fi. The details of this can be found in |7(. The scalar 
field is minimally coupled to Einstein gravity, non-minimally coupled to matter, and has a 
potential given by 



V{(j>) = p 2 Aff exp [ -2 J I 
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Consider vacuum cosmological solutions. We must specify the initial values of and <ft', 
denoted as <pi and 0^. For simplicity we take & -C Mp. There are three qualitatively distinct 
outcomes, depending on the value of 0' \. 

1. Eternal de Sitter. There is a critical value of <p' \ = <p' c for which cf) just reaches the 
maximum of the potential V(<ft) and comes to rest. In this case the Universe asymptotically 
evolves to a de Sitter solution. This solution requires tuning and is unstable, since any 
perturbation will induce the field to roll away from the maximum of its potential. 

2. Power-Law Acceleration. For (f)' \ > <p' c , the field overshoots the maximum of V(4>) 
and the Universe evolves to late-time power-law inflation, with observational consequences 
similar to dark energy with equation-of-state parameter w DE = —2/3. 

3. Future Singularity. For ' i < <p' c , <fi does not reach the maximum of its potential and 
rolls back down to <fi = 0. This yields a future curvature singularity. 

In the more interesting case in which the Universe contains matter, it is possible to show 
that the three possible cosmic futures identified in the vacuum case remain in the presence 
of matter. 

By choosing fi ~ 10~ 33 eV, the corrections to the standard cosmology only become im- 
portant at the present epoch, making this theory a candidate to explain the observed accel- 
eration of the Universe without recourse to dark energy. Since we have no particular reason 
for choosing this value of fi, such a tuning is certainly not attractive. However, it is worth 
commenting that this small correction to the action, unlike most small corrections in physics, 
is destined to be important as the Universe evolves. 



Clearly our choice of correction to the gravitational action can be generalized. Terms of 
the form -^ n+1 ^/R n , with n > 1, lead to similar late-time self acceleration, with behavior 
similar to a dark energy component with equation of state parameter 

2(n + 2) 

w cS = -1 + — — — — — — . 3 

3(2n + l)[n+ 1) 
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Therefore, such modifications can easily accommodate current observational bounds 
on the equation of state parameter —1.45 < wde < —0.74 (95% confidence level). In the 
asymptotic regime, n = 1 is ruled out at this level, while n > 2 is allowed; even n = 1 is 
permitted if we are near the top of the potential. 

In this paper we seek to extend this approach. For the actions considered in this paper, 
a tranformation between the matter and Einstein frames does not necessarily make sense. 
Therefore, we would like to analyze the dynamics in the matter frame itself. A technique for 
this analysis is presented here. 

As an example, we begin by applying this technique to the CDTT model. We will work 
in the matter frame with a spatially flat Robert son- Walker metric 

d s 2 = -dt 2 + a{tf dx 2 , (4) 

with a(t) being the scale factor. The time-time component of the field equations, the Fried- 
mann equation in this non-standard cosmology, is 

3H 2 J- ( 2HH + lhH 2 H + 2H 2 + 6H A ) = , ( 5 ) 

12(H + 2H 2 ) 3 V J M 2 ' K J 

where an overdot denotes differentiation with respect to cosmic time, and H = a/a. 

To perform a phase-space analysis of such equations we write if as a function of H. Let 

x = -H(t) 

y = H(t) , (6) 

so that 

£ dH dH dH dy 
dt dH dt dx 
In this way we may write (0), a third order equation in a(t), as a second order equation 

in H(t) and hence as a first order equation in y(x). The fact that the FRW equation for 



x 2 



any f(R) theory can be reduced to the first-order equation for the spatially flat case (second 
order in case of a non-zero spatial curvature) was first introduced in , 3(| . 

Many cosmologically interesting solutions, including accelerating ones, are power-law so- 
lutions of the form a(t) oc t p . In such cases H = —H 2 /p (i.e. y = —x 2 /p). In anticipation 
of finding such solutions as asymptotic solutions to our equations, we define a new function 
v(x) by 

v(x) - 

y 

with i; / 0. Power-law solutions in the asymptotic future are then easily identified as 
v(x) — > p = constant as H — > 0. 
Furthermore, if 

\H\ = \y\ -»• oo, then |u| -»• if x ^ 0. (9) 

So if x is not zero, as v — > we approach the singularity \H\ — > oo. This trick is invoked 
throughout this paper. 

Here, let us apply it to the simple case of (J5J). The relevant first order equation is 

x ^ = 2v + — A \x A (36 - 216w + 432t» 2 - 288w 3 ) + ^ (2v - 15v 2 + Qv 3 )] , (10) 
ax 2/x 4 

with the resulting phase plot shown in figure ^ 

Since the x-axis on this plot is (minus) the Hubble parameter, earlier times in the universe 
lie to the (negative) left and late times lie closer to x = (in all but exponential or phantom 
evolution). 

Note that the numerical solution shows the accelerating attractor a(t) oc t 2 , corresponding 
to v (x) — > 2 as x — > 0, as expected from the analytic, Einstein-frame method. Indeed, for 
general n this attractor, at v — (2n + l)(n + l)/(n + 2) can be obtained directly from the 
asymptotic form of the generalization of equation (jlOj) . In addition, our method pinpoints a 
singularity in the phase space, corresponding to a power-law evolution with exponent p = 1/2 
(that of radiation). Both features are also evident from the study of the asymptotically late 
time behavior of equation JTUJ). In order to have a constant v as a solution, we need 

2v 2 - 5wo + 2 = , (11) 
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FIG. 1: Two phase portraits for the modified gravity model proposed in jj. The left portrait is in 
the coordinates (x, v), for which an attractor at constant v = p corresponds to a power-law solution 
with a(t) oc t p . The right portrait is for the same theory in the (H,H) plane, with the unstable de 
Sitter solution at (0, 1). 



which has two real solutions: v = 1/2, 2, as expected. 

The singularity at p = 1/2 occurs because the Ricci scalar vanishes for this particular 
power. One might worry that this is problematic for describing the radiation-dominated 
phase in standard cosmology, since nucleosynthesis occurs during this epoch and provides a 
particularly strong constraint on deviations from the standard Friedmann equation at that 
time. However, as we shall see later, this is not a problem when matter sources are included 
explicitly, since even during radiation domination the Ricci scalar does not vanish exactly, 
but rather, has a small contribution from non-relativistic matter. 

Nevertheless, the singularity is a new feature that was not found in the Einstein-frame 
analysis Q. This is presumably because R = is a singular point of the conformal transfor- 
mation used to reach the Einstein frame. We shall see similar singularities for some of the 
more general actions we consider in this paper. 

Another way to visualize the solutions to our models is to use a more traditional phase 
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portrait in the (H, H) plane. 

The outline of this paper is as follows. In the next section we shall introduce the general 
class of actions we are interested in. In section ITTT1 we analyze the vacuum equations, describ- 
ing the singularity and attractor structure in detail before moving on to some simple special 
cases. In section I1VI we introduce matter into the equations, demonstrating briefly that the 
late-time attractor solutions of the system remain unchanged and setting up the formalism 
used in the appendix to establish stability of the system. In section IS we summarize our 
findings and comment on the status of these models as origins of cosmic acceleration. The 
paper contains two appendices. Appendix A contains definitions of a number of functions 
used in the body of the paper and Appendix B consists of a proof of the stability of the 
vacuum solutions under the addition of matter. 



II. A GENERAL NEW GRAVITATIONAL ACTION 

We now generalize the action of []| to include other curvature invariants. There are, 
of course, any number of terms that we could consider. We have chosen to consider those 
invariants of lowest mass dimension that are also parity-conserving 



P = R^R^ 

Q = Raf3~/8 R a ^ S 



(12) 



Since we are interested in adding terms to the action that explicitly forbid flat space as a 
solution, we will, in a similar way as in consider inverse powers of the above invariants. 

It is likely that such terms introduce ghost degrees of freedom. We shall not address this 
problem here, since it is beyond the scope of this paper. Rather, if ghosts arise we shall 
require that some as yet unknown mechanism (for example, extra-dimensional effects) cut 
off the theory in such a way that the associated instabilities do not appear on cosmological 
time scales js7| (see w| for an example of a concrete model where ghosts are brought under 
control by higher-derivative terms). We therefore consider actions of the form 



S= d 4 xy/=g[R + f(R,P,Q)]+ d^x^gC 



(13) 



where f(R, P, Q) is a general function describing deviations from general relativity. 
It is convenient to define 

f - 9f , _df _df 

fR= dR' fp= dP' fQ= dQ' (14) 
in terms of which the equations of motion are 

RfjLV 2 ^ 2 9fU> f 

+ fft RfJ,u + 2/p -R^ R av + 2/q Ra/3-yn R^ 1 v 

+ ^ D/h - V^V./p - 2V a Vp[f P R a { ^ u) ] + D(/ P R, v ) 

+ g^V a Vp{f P R al3 )-AV a Vp[f Q R a { ^} = SirGT^ . (15) 

It is straightforward to check that these equations reduce to those of the simple model 
of Q for f(R) = — /i 4 / R. 

We would like to obtain constant curvature vacuum solutions to these field equations. 
To do so, we take the trace of (|T5Jl and substitute Q = R 2 /4 and P = R 2 /6 (which are 
identities satisfied by constant curvature spacetimes) into the resulting equation to obtain 
the algebraic equation: 

(2/ Q + 3/ P ) R 2 + 6 (f R ~ 1) R ~ 12/ = . (16) 

Solving this equation for the Ricci scalar yields the constant curvature vacuum solutions. 

Evidently, actions of the form ()13|) generically admit a maximally-symmetric solution: 
R = a non-zero constant. However, an equally generic feature of such models is that this 
deSitter solution is unstable. In the CDTT model the instability is to an accelerating power- 
law attractor. This is a possibility that we will also see in many of the more general models 
under consideration here. 

Before we leave this section, note that, in analyzing the cosmology of these general models, 
it is useful to have at hand the following expressions, which hold in a flat FRW background 



R = 6(4 + -) =6(H + 2H 2 ) (17) 

P = 12(4 + 4 + ^-) =m(H + H 2 ) 2 + H* + H 2 (H + H 2 )\ (IS) 
ar a z a z a 

Q = 12(4 + 41 =12[(H + H 2 ) 2 + H 4 } . (19) 
or a z 



We have provided these both in terms of the scale factor a(t) and in terms of the Hubble 
parameter H(t) = a(t)/a(t), since they will be separately important in this paper. 



III. VACUUM SOLUTIONS 

In this section we study cosmological solutions to the field equations (fT5Jl in the absence of 
matter sources. Physically, this is important because we are hoping to find novel cosmological 
consequences arising purely from the gravitational sector of the theory. Mathematically, this 
provides us with valuable insight into the structure of the equations, which take a significantly 
simplified form wherein the Hubble parameter may be treated as the independent variable. 

As mentioned in the previous section, we will consider inverse powers of our curvature 
invariants and, for simplicity, we will specialize to a class of actions with 



where n is a positive integer, fi has dimensions of mass and a, b and c are dimensionless 
constants. M We will focus on the case n — 1 for most of the paper, because the analysis 
is less involved for that case. For general n the qualitative features of the system are as for 
n = 1 and we discuss the quantitative differences in our conclusions. 

A. Distinguished Points of the Action (and Equations) 

In general, the analogue of the Friedmann equation may be written in the following 
convenient form 



where A = A(H,H), B = B(H,H) and C = C(H,H) arise from the gravitational part of 
the action and M = M(a) describes the possible inclusion of matter. 

4 I} 

[1] Another potentially interesting possibility is a correction of form f(R,P,Q) = t-pr- This term has the 
same mass dimension as the /i 4 /-R term of CDTT. However, it turns out that this model does not possess 
any accelerating attractors. Specifically, the scale factor asymptotically approaches a(t) oc t 37 . 




4n+2 




(20) 
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It is also convenient to write this schematically in terms of our variables of the previous 
section as 

dv_ = x 6 f(v) +n 6 v 2 g(v) 
X dx 2^ v h(v) ' 1 ' 

where /(f), g(v) and h(v) are 6th, 4th and 2nd degree polynomials respectively in the variable 
v , whose explicit form is given in Appendix A. We are often interested in a particular subset 
of the phase space, the region where v > 1/2. This is because from nucleosynthesis to the 
matter-dominated epoch, we expect Einstein's equations to provide a good approximation 
to the dynamics, and therefore, when matter becomes subdominant we should have 1/2 < 
v < 2/3. 

There are three types of special points in the phase space plots of these equations: 

1. Singular Points of the Friedmann Equation 

In our introduction we reconsidered the model of [?J and discovered a singular point of 
both the action and the equations of motion, corresponding to a power-law evolution 
with exponent p = 1/2. In this particular case the singularity occurred because R = 
for p — 1/2. In our more general models, analogous singularities occur whenever the 
denominator of the Friedmann equation blows up; i.e. at zeros of C, where C is defined 
by f|2ip. For the flat cosmological ansatz this occurs when 

H + h2 = -t(^) , (23) 



where we have defined 



4 \ H 



12a + 46 + 4c , njI , 

a = ; . (24) 

12a + 36 + 2c v ; 



Recall that a, 6, and c are parameters in the Lagrangian (|20|1. In our variables of the 
previous section this becomes 

1 



V = Pli2 = -(l±v / r^) , (25) 



each zero having multiplicity 3. 
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These singularities only exist if px,2 are real, i.e. a < 1 and therefore there are many 
invariants that do not admit this type of singularity. In the simple case of corre- 
sponding to 6 = c = 0, note that we have a = 1 and recover p c = 1/2 as expected. If 
a < 1 

Pi < 1/2 p 2 > 1/2 (26) 

2. Singular Points at which H — > oo 

Points at which if — > oo occur at zeros of 5 and, in our variables of the previous 
section, correspond to 

+oo (27) 



dx 

at finite x and v. We denote the zeros by v(x) = v±,V2, where the Vi are (in general 
complex) constants constructed from a, b and c. 

If g(v) 7^ one singular point is at x = 0. Otherwise the singularities occur at solutions 
of v h(v) = 0, i.e. 

v[(108a 2 + hlab + 76c + 30ac + 2c 2 + 6b 2 ) v 2 

- (63a6 + 6c 2 + 9b 2 + 108a 2 + 156c + 54ac) v 

+ 18ab + 3c 2 + 66c + 27a 2 + 18ac + 36 2 ] = , (28) 

which are given by 



3a / la — 1 
v 3 = . (29) 



3. Late time Stable Points. 

Finally, we look for late time power-law attractors. This happens if v(x) —>■ constant 
(distinct from our previously mentioned singular values). Taking an asymptotic limit 
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of the equations of motion yields 

v [(288a 2 + 186 2 + 8c 2 + lUab + 96ac + 246c) v 4 

- (1512a 2 + 906 2 + 36c 2 + 738a6 + 468ac + 1146c) v 3 

+ (1836a 2 + 1236 2 + 62c 2 + 951a6 + 678ac + 1756c) v 2 
-(846a 2 + 696 2 + 44c 2 + 489a6 + 414ac + 1136c) v 

+ 135a 2 + 156 2 + 15c 2 + 90a6 + 90ac + 306c] = , (30) 

with solutions 



20 - 3a / / 120a 
si,2 = o U± a/1 



8 V y (20 -3a) 2 
S3 = Pi (31) 
s 4 = pi (32) 
s 5 = , (33) 

and s 2 > Si- 
It is clear that only s 1 and s 2 can be late-time power-law solutions. The other ones 
represent singularities; v — implies that y = H — > oo, whereas s 3 and s 4 are the 
singular points we discussed earlier. 

B. Summary of Possibilities 

Here we summarize the different vacuum possibilities. It is useful to define the following 
two constants 

280 + 60^ „ , 4X 

m = 1= « 1.78 34 

111 + 60^ 

280 - 60^ , . 

V2 = ^= « 24.9 . 35 

111 - 60V3 

1. a < 1 . In this case i>j are complex, whereas Pi and Sj are real. It is straightforward to 
show that s 2 > £>2 > 1/2 and s 2 > 1. 
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• < a < 1. In this case 1/2 > s± > p\ > and 1 > p 2 > 1/2- Solutions close to 
P2 are repelled from it, whereas si is an attractor. This leads to decelerating late 
time behavior. 

• 0>a>— 160/9. Here s 1 < pi < and solutions are attracted to (x — 0,v — 0). 
Furthermore p 2 > 1 . 

• a < — 160/9. In this case pi < si < and again we have that solutions are 
attracted to (x = 0,v = 0). Again, p 2 > 1- 

2. a > 1 . In this case V; L are real, whereas p; L are complex. 

• 1 < a < 4/3. Here vi < 1/2 and 1/2 < si < v 2 < 1. Furthermore s 2 > 1. Both 
Si and s 2 ar e attractors, but si describes a decelerating phase. 

• 4/3 < a < r\\. Here 1/2 < f 1 < 1 and f 2 > 1. Both (x = 0,v = 0) and s 1)2 are 
attractors, with the separatrix being at t>i j2 . 

• rji < a < 4. In this case 1/2 < v\ < 1 and f 2 > 1. Again, both (x = 0, v — 0) 
and V2 are attractors and the separatrix is at v±. There are no real solutions for 

• 4 < a < r/2. Here there are no real late-time attractors since the Sj are complex. 
We also have v 2 < and v± > 1. Solutions either evolve to (x = 0,t> = 0) or to 
v = +oo, with separatrix at v±. 

• a > r] 2 . In this case t> 2 < 0, v i > 1 and si < s 2 < 0. Evolution is to a decelerating 
attractor. 

3. a = 1 . This yields a promising class of solutions. We have 



never encounter this point. It is simple to show that solutions evolve to a late-time 
power-law attractor describing an accelerating phase, with 



1 

Pi = P2 = vi = v 2 = Si = - . 



(36) 



All singularities occur as v(x) 



1/2. However, from nucleosynthesis onwards we 



s 2 = —r 




(37) 
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which is otherwise independent of a, 6, c. 

The results of this section may be summarized in figure showing the values of the 
various distinguished points as a is varied. 




FIG. 2: The values of the various distinguished points as a is varied. 



C. Inverse Powers of P = R^ v R^ v 

Let us begin by dealing only with actions containing modifications involving P = R fiV . 
Our prototype is to consider f(P) = —m 6 /P, with m a parameter with dimensions of mass. 

Using (fTo^) we can see that there is a constant curvature vacuum solution to this action 
given by 

i£L = (16) 1/3 m 2 . (38) 

However, we would like to investigate other cosmological solutions and analyze their stability. 

From p5j1 . with the flat cosmological ansatz, the analogue of the Friedmann equation 
becomes 

m 6 



3H Z - 



8{3H* + 3H 2 H + H 2 f 



+ 33# 4 # 2 + 30H 6 H + 6H 3 HH + 6H 8 + 4H 5 H 



H A + 11H 2 H 3 + 2HH 2 H 

= . (39) 
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We analyze this equation using the same technique, with the same definitions, as in the 
example of the previous section. The relevant equation is 

dV _ n 1 



x 



dx 



2v 2 + 



- [-x 6 (24 



216v + 8QAv A - 19Uv a 



2m 6 v(2v 2 -3v + 

+ 2592v 4 - 19Uv 5 + 648v e ) + m 6 (v 2 - llv 3 + 33v 4 - 30v 5 + 6v 6 )] . (40) 



The solution to this equation is displayed graphically in figureEl We identify four fixed points 
of the system; two attractors at v ~ 0.77 and v ~ 3.22 and two repellers at v ~ 0.5 and 
v = 1. Clearly, in order to obtain a late-time accelerating solution (p > 1), it is necessary to 
give accelerating initial conditions (d > 0), otherwise the system is in the basin of attraction 
of the non-accelerating attractor at p ~ 0.77. 
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FIG. 3: Two phase portraits for the f(R,P,Q) = —m 6 /P modification. The left portrait is in the 
coordinates (x,v), for which an attractor at constant v = p corresponds to a power-law solution 
with a(t) oc t p . The right portrait is for the same theory in the (H, H) plane. There are two 
late-time power law attractors corresponding to p = (4 ± \/6)/2. The (-) branch, represented by 
the solid lines, is non-accelerating, while the (+) branch, represented by the dash-dotted line, is 
accelerating. 
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The exact exponents of the two late-time attractors of the system are obtained by studying 
the asymptotic behavior of ([40)1 . Substituting in a power-law ansatz and taking the late-time 
limit we find that, in order to have a constant v = v o as a solution, the exponent must satisfy 

Gvq - 30ug + 41v% - 23v + 5 = 0. (41) 

This equation has two real solutions (and two complex ones). The real solutions are: vq = 
2 - y/G/2 ~ 0.77 and v Q = 2 + y/E/2 ~ 3.22. 

Even the non-accelerating attractor is of some interest in this model. In order for structure 
to form in the universe, there must be a sufficiently long epoch of matter domination, for 
which a(t) oc t 2 ^ 3 . As matter redshifts away, however, since the universe is decelerating, we 
expect the universe to approach the attractor at p ~ 0.77. This corresponds to an effective 
equation of state w e R ~ —0.13; i.e. negative pressure, although not negative enough to 
provide a good fit to the supernova observations. 



D. Inverse Powers of Q = R a ^ 5 R 01 ^ 5 

Now let us move on to actions containing modifications involving only Q = R a /3-ys R a,3 ^ S - 
Our prototype example is f(Q) = —M 6 /Q, with M another parameter with dimensions 
of mass, and the analysis follows much the same as in the previous subsection, albeit with 
different results. 

Again, 1)16)1 demonstrates that there is a constant curvature vacuum solution to this action 
given by 

i?S st = (24) 1 / 3 M 2 . (42) 

What about other possible cosmological solutions? From (J15j) . with the flat cosmological 
ansatz, the analogue of the Friedmann equation becomes 



24[H 2 + 2H 2 H + 2H 4 } 3 

+ 20H 2 H 3 + 3H A + AH 5 H + 12H 3 HH + 6HH 2 H 



8H S + 3QH 6 H + 5AH 4 H 2 

. (43) 
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Employing our phase space technique once more, in the variables best-suited for analyzing 
power-law behavior, the relevant equation is 

x ^ = 2v + —- r } [-s 6 (576i; 6 - 1728v 5 + 2592v 4 - 2304v 3 

dx 2M G v[2v 2 - Qv + 3] L v 

+ 1296w 2 - 432w + 72) + M 6 (8v 6 - 36v 5 + 54v 4 - 20v 3 + 3v 2 )] . (44) 

It is clear from this that our action consisting of f(R, P, Q) = —M 6 /Q does not admit any 
late-time power-law attractors. 
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FIG. 4: Phase portrait for the f(R,P,Q) = —M^/Q modification in the coordinates [x, v), for 
which an attractor at constant v = p corresponds to a power-law solution with a(t) oc t p . 

This is consistent with a study of the late-time asymptotics of ([44)1 . As in our previous 
analysis, late-time power-law solutions correspond to real solutions to the equation 

8w 4 - 36v% + 62v% - Uv + 15 = . (45) 

However, no such real solutions exist, confirming our phase-space analysis. 
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IV. INCLUDING MATTER 



We now show that the late-time behavior of our vacuum solutions remains unaltered upon 
the inclusion of matter. 

We begin by rewriting the equation of motion in a more convenient form. Let be the 
tensor defined by the left hand side of (|15|). Then, the generalized Friedmann equation takes 
the form 

S 00 = 8ttG P , (46) 
where p is the energy density of a perfect fluid with equation of state p = wp. Now, since 



x 



-H, y = —x, and p oc a 3 ( 1+u, ) ? we have = ^(1 + w)-. Combining this relation 

with (}4l)j) yields the equation of motion 

x^ = 3(l + w)vZ m . (47) 
ax 

Thus far, we have not imposed any dynamics. Specializing to a theory with f(R, P, Q) = 
p 6 /(aR 2 + bP + cQ) gives 

Z 00 = ^F(x,v,s), (48) 

with 

where s — The explicit forms of the functions h(v),g(v), and d(v), defined in Appendix 
A, are not needed in this section. Substituting (|4*K|) and (|4H|) into (|4Tjl gives the equation 
governing the dynamics of our theory in the presence of matter: 

dF 

x—= 1 {v)F, (50) 
where 7(f) = [3(1 + w) v + 4] . Using the chain rule this becomes 

x{F x + F v s + F s s') = 7 (v) F , (51) 
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where 

F - = 1Z = *««"-V^ (52) 

<9F 6 v 2 g v + 2vg — 2xhs — 2xvh v s 6 v 2 g — 2xvhs 

v = ~fo =fl & ^ " 

<9F 2fi 6 xvh 

and a prime denotes differentiation with respect to x. 

Seeking late-time power law behavior, we take the limits x — > and s — > in ()51|) . 
This yields the condition g(t>o) = 0, which is solved by the same power law fixed points as 
those obtained in vacuum. This is as one might expect, since matter redshifts away in the 
asymptotic future. However, the above description proves useful for dealing with the issue of 
stability. Since this is somewhat technical, we relegate the details to Appendix B and merely 
assert here that these fixed points do remain stable in the presence of matter sources. 



V. COMMENTS AND CONCLUSIONS 



In the extreme low curvature regime, our only tests of general relativity are cosmological. 
The discovery of new phenomena at these scales may point to new matter sources, but 
alternatively may hint at hitherto undetected modifications of gravity. 

The acceleration of the universe provides a particular challenge to modifications of gravity. 
Unlike the known perturbative corrections to the Einstein-Hilbert action arising from string 
theory, late time acceleration requires modifications that become important at extremely 
low energies, so low that only today, at the largest scales in the universe, is the resulting 
curvature low enough to lead to measurable deviations from general relativity. 

These considerations led some of the current authors and others to consider new terms in 
the action for gravity that consist of inverse powers of the Ricci scalar. It is easy to show that 
such an approach introduces a de Sitter solution. However, this solution is unstable to a late- 
time accelerating power law attractor. For appropriate choices of parameters, this theory 
is a candidate to explain cosmic acceleration without the need for dark energy, although 
the simplest such theories are in conflict with solar system tests. For Lagrangians that are 
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functions of only the Ricci scalar, there exists a map to an Einstein frame, in which the 
new degrees of freedom are represented by a scalar field. As a result, such modified gravity 
theories share many features in common with some dark energy models. 

In this paper we have introduced a much more general class of modifications to the 
Einstein-Hilbert action, becoming relevant at extremely low curvatures. Specifically, we 
have considered inverse powers of arbitrary linear combinations of the curvature invariants 
R 2 , P = R [W R IXV and Q = R a /3yS R a/3lS . Such modifications are not simply equivalent to 
Einstein gravity plus scalar matter sources. 

We have performed a general analysis of the late-time evolution of cosmological solutions 
to these theories. Many of the theories exhibit late-time attractors of the form a(t) oc t p , 
with p some constant power. Indeed, there are often multiple such attractors. For a large 
class of theories there exists at least one attractor satisfying p > 1, corresponding to cosmic 
acceleration. 

The detailed structure of cosmological solutions to these theories turns out to be quite 
rich and varied, depending on the dimensionless parameters entering the particular linear 
combination. Two distinct types of singularities may exist, as well as the late-time power 
law attractors. We have identified all those theories for which the late time behavior is 
consistent with the observed acceleration of the universe, providing a whole new class of 
theories - generalized modified gravity theories - which are alternatives to dark energy. 

The results we have found for the modification / = fi 6 /(aR 2 + bP + cQ) may be gener- 
alized to modifications 



using exactly the same calculational techniques. 

In general there are power-law attractors with the following exponents, whenever they 



./ 



(aR 2 + bP + cQ) 



n 



(54) 



are real 




8n 2 + Wn + 2 - 3a ± y/T 
4(n + 1) 



(55) 



where 



r = 9n 2 a 2 - (80n 3 + 116n 2 + 40n + 4) a 

+ 64n 4 + 160n 3 + 132n 2 + 40n + 4 . 



(56) 
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Clearly, as n — > oo, the smaller attractor tends to 0, whereas the larger one increases linearly 
as An. 



Two special cases are 

/ = ' 
/ = ' 

for which the power law attractors, V\ 2 are 



m 4n+2 



pn ' 
M 4n+2 

(57) 



(P) I2n 2 + 9n + 3 ± ^UAn 2 + 120n 3 - 15n 2 - 30n - 3 
;i ' 2 ~ 2(3 + 3n) 



(5? 



f 0) 4n 2 + 2n + 1 ± Vl6n 4 - 16n 2 - lOn - 1 , . 

v\V = . (59) 

1,2 2n + 2 v ; 

For n — 1, all the above expressions agree with the values found earlier. It is interesting 
to note that are imaginary for n = 1, but are real for all n > 1. 

Of course, much remains to be done. We have not addressed solar system tests of these 
theories since this is a complicated analysis that is beyond the scope of our current paper. 
We have not focused on detailed comparisons between our models and the supernova data 
and it is possible that there are specific signatures of this new physics in such data. For 
example, the modified Friedmann equations arising from the theories presented here should 
directly provide information about the jerk parameter. We intend to consider such effects, 
not only in the present models, but in those proposed by us and other authors, in a separate 
paper focused on the connections of these models with observations. In this context, as with 
the case of the simple modifications introduced in []]], it may also be interesting to study 
more complicated functions of the curvature invariants we have considered. 
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APPENDIX A: SOME DEFINITIONS 



In Section mil we defined 



4(3a + 6 + c) 

a = , (Al) 

where 

A = 12a + 36 + 2c . (A2) 
We should distinguish betweeen two cases: 



1. A = 

In this case a diverges, but we may define 

g{v) = (3a + 6 + c) (282a + 696 + 44c) v + 15 (3a + 6 + cf (A3) 
h(v) = 3(3a + 6 + c) 2 (A4) 
d( v ) = 3a + 6 + c . (A5) 

Note that both A, and 3a+6+c, cannnot vanish simultaneously, as this causes the Friedmann 
equation to be singular for all values of x and v. 



2. A / 

In this case we have 



g(v) = 2A 2 (v- s 2 ) (v - Sl ) (v - p 2 ) (v - Pl ) (A6) 
h(v) = A 2 (l - (v - v 2 ) (v - v x ) (A7) 
d(v) = A(v- P2 ) {v-pi) , (A8) 
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where either f 12 or pi )2 are complex unless a = 1, in which case, = Pi,2 = si = 1/2 and 
S2 = 3.75. Also Si 5 2 in general may be complex. As we saw in Section ITTT1 Sj, Vi and pi are 
all functions of a only. 



APPENDIX B: STABILITY OF FIXED POINTS IN THE PRESENCE OF MAT- 
TER 



Let us rewrite ()51|) as a system of first-order ordinary differential equations 



x 



j(v) F F v 



xF s 
s . 



(Bl) 
(B2) 
(B3) 



Treating x as a dependent variable makes the system autonomous, facilitating a phase-space 
analysis. 

To obtain the phase portrait of the system we define a vector field by W T = (x', s', v') 
and plot this in the vicinity of the fixed point. In figure El we show a 2D slice of this phase 
portrait by choosing a section at constant x. Note that at the fixed point, = (1, 0, 0). 

To analyze the stability of the system, we first linearize the system about the fixed point 
v = v , s = 0, and x = x < 1 (all expressions evaluated at these values of v, s and x will 
carry the superscript ^) 

We write the linearized system of equations as 

a 

the matrix M being defined by 
and 



w. 



(o) 



(B4) 



1(0) 
J m 



1(0) 



(B5) 



r/ 



s 



(B6) 
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FIG. 5: Phase plot for the linearized equations close to the power law solution in the coordinates 
(x, v), for which an attractor at constant v = p corresponds to a power-law solution with a(t) oc t p . 

Now, introducing equilibrium coordinates 

Ml = X — Xq (B7) 
u 2 = s (B8) 
m 3 = V-Vq , (B9) 

we finally obtain the linearized equations 

u[ = 1 (BIO) 

, 7o + o- -l 70CT0 mn \ 

U 2 = U 2 2~M 3 (Bll) 

Xq Xq 

u' 3 = u 2 , (B12) 



With (T - V09v '° 



2 h 
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Notice that the first equation has decoupled from the other two. Hence, it suffices to 
study the behavior of the subsystem (u2,uz). We will use standard results from the theory 
of dynamical systems to do so (see 38]). 

In the vicinity of the fixed point, the behavior of the system can be classified by the 
eigenvalues of the sub-matrix M^-, with i,j = 2,3, with characteristic equation 



x 



M. 



(0) 



2v 



. 



(B13) 



Since we are only interested in values of v and w (the equation of state parameter) which 

nrr 



satisfy v > 
choice, and 



J2 and — 1 < w < 1 (see laa |4d, l4l| for arguments why this is a sensible 
for how one may be tricked into inferring values outside this range 



when considering gravitational theories other than General Relativity), we have that 7 > 1. 
Furthermore if we choose Vq = S2 (see Section HIT]) , the value of the bigger power law, typically 
o"o > 0. This implies that both M^J and M^j are negative, so that s 2 is either a stable node 
or a stable spiraling helix. The latter case is tantamount to stability. 
We have a stable spiralizing helix if 



(M 



(Oh 2 



+ 4 M 2 ( ? < 



If not, we get a stable node. Relation (|B14|) leads to the condition 

(7o + cr - I) 2 - 47 cr < , 



(B14) 



(B15) 



which holds for all Xq <C 1. 

Applying the last condition to the case a = 1 and w = (dust), we have that (see 
Appendix |XJ) 

1 (B16) 



= 3^2(52 -0.5) = 16.25 



and 



7o = 3s 2 + 4 = 15.25. (B17) 
For these values of a and 7, (IB 15 J) is easily satisfied, so that the power law solution s 2 = 3.75 
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is a stable spiraling helix. 
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